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In order to describe heavy-ion fusion reactions around the Coulomb barrier with an actinide tar- 
get nucleus, we propose a model which combines the coupled-channels approach and a fluctuation- 
dissipation model for dynamical calculations. This model takes into account couplings to the col- 
lective states of the interacting nuclei in the penetration of the Coulomb barrier and the subsequent 
dynamical evolution of a nuclear shape from the contact configuration. In the fluctuation-dissipation 
model with a Langevin equation, the effect of nuclear orientation at the initial impact on the pro- 
lately deformed target nucleus is considered. Fusion-flssion, quasi-flssion and deep quasi-flssion are 
separated as different Langevin trajectories on the potential energy surface. Using this model, we 
analyze the experimental data for the mass distribution of flssion fragments (MDFF) in the reac- 
tions of 3*'36g^238ij ^^^ 30g^_|_238][j ^^ several incident energies around the Coulomb barrier. We 
flnd that the time scale in the quasi-flssion as well as the deformation of flssion fragments at the 
scission point are different between the 30gj_|_238 
asymmetries of the quasi-flssion. 



U and 36s_|-238u systems, causing different mass 
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I. INTRODUCTION 

The prediction of the existence of the "Island of Stabil- 
ity" in the nuclear chart has encouraged searches of new 
elements [i| . The synthesis of these superheavy elements 
has been carried out using heavy-ion fusion reactions be- 
tween stable nuclei, in which two different types of reac- 
tion have been employed. In the cold fusion reactions, 
lead and bismuth targets are used [2, Q- The element 
with Z = 113 reported by RIKEN used this type of re- 
action Q. The superheavy nuclei (SHN) synthesized in 
the cold fusion reaction produce nuclei with relatively 
small number of neutrons. The other type of reaction, 
called the hot fusion reaction, on the other hand, uses 
actinide nuclei as targets. With this type of reaction, 
production of elements with Z — 114, 115, 116, 117 and 
118 were reported by the Flerov Laboratory of Nuclear 
Reactions (FLNR) [Sj . These nuclei as well as those pro- 
duced as descendants in the a-decay chain have relatively 
larger number of neutrons. Recently, other laboratories 
than FLNR also performed experiments of hot fusion re- 
actions and obtained results that are consistent with the 
data by FLNR p-Q. At present, attempts to produce 
elements 119 and 120 are made or planed in several fa- 
cilities using actinide targets. In order to produce new 
nuclei or elements not discovered so far, an accurate pre- 
diction of the production cross sections is an important 
issues in the SHN research. 

An important quantity for a prediction of the cross 
section is the probability of fusion after the interacting 
nuclei have the initial contact. Due to the complexity of 
the process, however, a good method to predict the fusion 



probability has not been established well. Furthermore, 
actinide nuclei are prolately deformed, so that one needs 
to introduce nuclear orientation as an additional degree 
of freedom. The effect of nuclear orientation on fusion 
probability has been evident already in the reactions us- 
ing rare-earth nuclei with a prolate shape [iCt lllj . This 
makes the hot fusion in contrast to the cold fusion reac- 
tions, which use spherical nuclei. 

Our strategy to calculate the fusion probability is to 
use the unified model I13|. which has been developed by 
the FLNR theory group |13 - [l4| . The model can describe 
every entrance and exit channels in heavy-ion collisions, 
and can calculate also the time evolution of the nuclear 
shape, where the binary decay of the composite system 
(fission) can be treated. In this model, a system first 
feels a diabatic potential in the early stage of the initial 
collision. The potential is then gradually shifted to an 
adiabatic potential in a time-dependent manner. A tra- 
jectory calculation is performed on the time-dependent 
unified potential energy surface by using the Langevin 
equation. In the trajectory analysis, different types of 
fission can be separated, that is, fusion-fission (FF) and 
quasi-fission process (QF). The fusion- fission is a fission 
of a compound nucleus, and is defined as the case when 
the trajectory enters the region of compound-nucleus, fol- 
lowed by fission. The quasi-fission, on the other hand, is 
the fission event whose trajectory does not enter the re- 
gion of compound nucleus. The fusion probability is de- 
fined as the FF events normalized to all the fission events 
(that is, FF-I-QF). The model should be constrained or 
checked by experiments by investigating whether the cal- 
culation can reproduce the measured fission spectra such 



as mass and total kinetic energy distributions. 

Recently, the mass distributions of fission fragments 
(MDFF) for the reactions 36'34g^238u ^^^^ 30gj_,_238u 
at several incident energies around the Coulomb bar- 
rier were measured by the Japan Atomic Energy Agency 
(JAEA) |l5l4l8l |. One of the findings in the experiment 
is that the mass asymmetry in QF is different between 
the 3°Si+^^*U and 36.34g_^238u systems at low incident 
energies. In this work, we attempt to analyze these data 
in order to understand the reaction mechanism. 

In the previous paper [l^ , to estimate the capture and 
fusion cross sections, we only considered the spherical- 
spherical configuration as the first approximation, limit- 
ing to the energy region above the Bass barrier. In or- 
der to extend the calculation down to energies below the 
barrier, the effect of nuclear orientation has to be taken 
into account. Such effect has been well established in 
the approaching phase of the reaction using the coupled- 
channels approach |2(]| - |22| . However, it is still difficult 
to calculate the adiabatic potential energy to be used in 
the unified model with the two-center parametrization for 
subsequent shapes of the nuclear system, starting from 
the configuration of arbitrarily oriented two deformed nu- 
clei touching each other to the spherical compound nu- 
clei. In this paper, we propose a new model which can 
avoid this difficulty. In the new model, all the orientation 
angles are effectively taken into account by introducing 
the effective charge-center distance at the contact point 
as a variable. With this prescription, a dynamical cal- 
culation for superheavy elements is possible for the first 
time at energies below the Coulomb barrier. 

The paper is organized as follows. In Sec. II, we detail 
the framework of the new model, which combines the 



coupled-channels method and the dynamical Langevin 
calculation. In Sec. Ill, we show the results for the cross 
sections and MDFF at several incident energies for the 



reactions of 



36,34t 



236 



U and 30Si 



236 



U. In Sec. IV, 



we discuss the reasoning of the different shapes of MDFF 
observed in these reactions. In Sec. V, we present a sum- 
mary of this study and further discussion. 



II. MODEL 

A. Coupled-channels method 

Excitations of the rotational states in a deformed nu- 
cleus in the approaching phase of heavy-ion collisions 
considerably modify the fusion barrier and increase the 
capture cross sections at the sub-barrier energies. For 
a heavy deformed nucleus considered in this paper, it is 
reasonable to introduce the sudden approximation to the 
coupled-channels (CC) equations. In this approximation, 
the capture cross section is given as |20l - |22| . 
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where 6 is the angle of the incident projectile nucleus with 
respect to the symmetry axis of the deformed target, and 
Ecm denotes the incident energy in the center-of-mass 
frame. Here, we have assumed that the target nucleus has 
an axially symmetric shape. (Tcap(-Ecm; ^) is the capture 
cross section for a given value of 0, calculated with an 
angle dependent internucleus potential Vcc(^j^)j 
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Here, Rt is the equivalent sharp surface radius of the tar- 
get nucleus, and P2 and ^4 stand for the quadrupole and 
hexadecapole deformation parameters of the target nu- 
cleus, respectively. We have assumed the Woods-Saxon 
function for the nuclear potential, V^^K Using the pene- 
tration probability for the i-th partial wave, Ti, the cap- 
ture cross sections are given by 



I 

taken into account for each 9 in the same way as in the 
computer code CCFULL (23| . 

The fusion cross section is calculated by multiplying 
the probability to form a compound nucleus, Pcn, to the 
capture probability, Ti(Ecm',9), at each incident angle 9 
and integrating it over the solid angle as 
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where k is the wave number of the incident flux. The 
vibrational excitations in the projectile nucleus are also 



CTfusiEcm) = / rf(cOs6')crfus(£'cm;6'), 
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with 
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As we explain in the next subsection, the formation prob- 
abihty Pqn is estimated with the dynamical calculation 
using the Langevin equation. In order to introduce the 
orientation angle dependence to Pcn{E,£,9), we make 
the following approximation. We first notice that the 
charge-center distance at the touching point is strongly 
correlated with the incident angle 9. In our model, the 
touching distance is assumed to be 

^touch(^) = R + RtP2Y2o{0) + RrPiYioie), (8) 

where R is the radius parameter in the Woods-Saxon po- 
tential, Eq. ^. In the dynamical calculation, we use an 
angle independent potential energy surface, but starting 
from the angle dependent initial touching distance given 
by Eq. ([8]). With this prescription, we can for the first 
time extend the dynamical Langevin calculation down to 
the subbarrier region. 



B. Dynamical calculation 

After the projectile enters with an arbitrary orienta- 
tion relative to the symmetry axis of the deformed tar- 
get nucleus, the collision is replaced by the one from the 
polar-side of the target nucleus and the trajectory calcu- 
lation starts from the configuration corresponding to the 
touching distance ztouch- That is, we consider only the 
nose-to-nose configuration. In this stage, we assume that 
the potential has been shifted to the adiabatic potential 
from the diabatic one. In the reactions of ^'^^U, the static 
deformation of /Sa == 0.215 (6 - 0.2) [24J is used. 

The nuclear shape is defined by the two-center 
parametrization [2^, [2Q|, which has three deformation 
parameters, zo,S, and a. zq is the distance between two 
potential centers, while a = {Ai — A2)/{Ai + A2) is the 
mass asymmetry of the colliding nuclei, where Ai and 
A2 denote the mass numbers of heavy and light nuclei, 
respectively [27| . 6 denotes the deformation of the frag- 
ments, and is defined hy S — 3{R\\ — Ri_)/{2R\\ + -Rj_), 
where R\\ and R± are the half length of the axes of an 
ellipse in the zq and p directions of the cylindrical coor- 
dinate, respectively, as shown in Fig. 1 in Ref. 



251. We 



assume that each fragment has the same deformation. 
The deformation parameters 5 and /32 are related to each 
other as 



ti2 = 



iki^~s) 



(9) 



Notice S < 1.5, because i?j| > and R± > 0. In order to 
reduce the computational time, we employ the coordinate 
z defined as z = zq/{RcnB), where Rcn denotes the 



radius of a spherical compound nucleus and B is defined 
asS = (3-H(5)/(3-2(5). 

The neck parameter e entering in the two-center 
parametrization has been adjusted in Ref. [2^ to re- 
produce the available data, assuming different values be- 
tween the entrance and the exit channels of the reactions. 
In our study, we use e = 1 for the entrance channel and 
e = 0.35 for the exit channel. We assume the follow- 
ing time dependence for the adiabatic potential [23], ex- 
pressed in terms of the relaxation time t^ for e: 

V{q, t) - V{q, e ^ l)/,(t) + V{q, e - 0.35)[1 - /,(i)], 

(10) 
with 
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Here, q — {z, d, a} is the deformation coordinate. We use 
Te — 10^^° sec, as it gives a reasonable account for the 
available data of mass, angular, and kinetic energy dis- 
tributions of outgoing fragments produced in heavy-ion 
collisions (that is, deep-inelastic reaction, nuclear trans- 
fer reactions, fusion- fission) |12h141 [19|. 

For a given value of e and a temperature of a system, 
T, the adiabatic potential energy is defined as 
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Here, Vld is the potential energy calculated with the 
finite-range liquid drop model, given as a sum of of the 
surface energy Es 3Q\ and the Coulomb energy Eq ■ Vsh 
is the shell correction energy evaluated for each tempera- 
ture using the factor <I>(T), in which E^ is the shell damp- 
ing energy chosen to be 20 MeV [3l| and a is the level 
density parameter. At the zero temperature (T = 0), 
the shell correction energy reduces to -E^hoii- ^^^ second 
term on the right hand side of Eq. (|12l) is the rotational 
energy for an angular momentum £ [27| . with a moment 
of inertia, I{q)- 

Since we employ the CC model for the approaching 
phase to describe the penetration of the Coulomb barrier, 
the two-body part of the unified mode is omitted |13 - [l4l 
Il9| in our calculations. The multidimensional Langevin 
equations [l^, [23| are thus simplified as 
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where pi = dqi/dt with i = {z, 5, a}. The summation is 
performed over repeated indices. In the Lengevin equa- 
tion, ra-ij and 7ij are the shape-dependent coUective iner- 
tia parameter and the friction tensor, respectively. The 
wall- and- window one-body dissipation [32r;34j is adopted 
for the friction tensor. A hydro dynamical inertia ten- 
sor is adopted with the Werner- Wheeler approximation 
for the velocity field [35j. The normalized random force 
Ri{t) is assumed to be white noise, i.e., {i?j(i))=0 and 
{Ri{ti)Rj{t2)) = 25ij5{ti - ia). The strength of the ran- 
dom force Qij is given by 7yT = Y.k 9ij9jk- 

The temperature T is calculated from the intrinsic en- 
ergy of the composite system as Ei^t — aT^, where i^int 
is calculated at each step of a trajectory calculation as 
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The excitation energy of the compound nucleus E* is 
given by E* = E^ra — Q, where Q denotes the Q- value of 
the reaction. 

The fusion probability Pqn in Eq. (7) is determined 
in our model calculation by identifying the different tra- 
jectories on the deformation space. It is equivalent to 
the number of trajectories of compound-nucleus fission 
normalized to all the fission events. Formation of the 
compound nucleus is defined as the case that a trajec- 
tory enters in a compact-shape region in the adiabatic 
potential energy surface. We define the compound nu- 
cleus region (that is, the fusion box) by referring to the 
ridge of the fission barrier in the coordinate space |27'l . 
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FIG. 1. (Color online) Excitation functions of (Tcap, crA/2±20 
and cTfus for the 36s_)_238u j-g^g^JQjj 'pjjg experimental data 
of (Tfiss, denoted by the circles, are taken from Ref. [isl ]. 



III. MASS DISTRIBUTION OF FISSION 
FRAGMENTS AND CROSS SECTIONS 



Reaction of ' 



lg_^238u 



Recently MDFF and the fission cross sections (crfiss) 
for the the 36,34g _|_ 238^ reactions were measured by the 
JAEA group [la, Il3]- In the experiment, fission events 
were selected in which the momentum of projectile is fully 
transferred to the composite system. For these systems, 
the fission cross sections are almost equal to those of the 
projectiles captured inside the Coulomb barrier, (Tcap- 

We first analyze these systems with the new model 
proposed in the previous section. The dashed-dot curve 
in Fig. [1] shows the calculated capture cross sections 
based on the coupled-channels model for the the reaction 
36g _j_ 238jj g^g g^ function of the incident energy. The 
experimental data are taken from Ref. jl5|. As in Refs. 
[11 [H-iil, we use the computer code CCDEGEN ;39] by 
taking into account the static deformation for ^38^ -^^jth 
(/32, /34)= (0.275,0.05). Couplings to the 2+ vibrational 
state at 3.29 MeV in ^^S (with ^2 =0.61 [13) and the 
3" state at 0.73 MeV in ^^^V (with /?3 = 0.086 [4l|) are 
also considered. We use ^0=105.0 MeV, i?=10.92 fm, 
i?T=7. 44 fm, and aws =0.75 fm for the Woods-Saxon 
potential, V^-'^''{r,9), given by Eq. ([3]). One can see that 



this calculation reproduces the measured cross sections 
down to the lowest incident energy below the Bass barrier 
(FBa.s = 158.8 MeV) [H. 

The solid curve in Fig. [T] shows the fusion cross sec- 
tion CTfus obtained by the new model with CC and the 
Langevin equation. At the Bass barrier, we obtain 
^CN = 0.03 for i = and 6 = 0. The dashed line shown 
in Fig. [T] denotes the cross section cr^/2±20: which is de- 
rived from the yield of the fission fragments whose mass 
number is located within ±20 around the symmetric fis- 
sion Acn/2. Notice that the fusion cross sections CTfus are 
significantly smaller than <yA/2±20- This indicates that 
the mass symmetric fission does not necessarily originate 
from the compound-nucleus state. 

In the previous study [l9|, we started the dynamical 
calculation with the spherical-spherical configuration in 
order to estimate the fission cross sections of the various 
kinds. Because the calculation using the Langevin equa- 
tion is a classical one, this did not allow us to calculate 
the cross sections below the barrier. With the new ap- 
proach, by considering the nuclear shapes at the contact 
configuration for each orientation, we can now obtain the 
cross sections also below the Bass barrier region. 

Above the Bass barrier region, the calculated CTfus with 
the present model is about 7 times larger than that in the 
previous study |l9|. In the present model, the Langevin 
calculation is started at the touching point assuming all 
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FIG. 2. (Color online) Mass distributions of fission fragments 
for the reaction of ^^S+'^^*U. The experimental data and the 
calculated results are denoted by the circles [13 and the his- 
tograms, respectively. The shaded areas show the calculated 
fusion-fission events. 



The results for the MDFF for the reaction of ^^S 
+ 238u g^j.g compared with the experimental data [l^ 
in Fig. [2] at seven incident energies from _Ecm= 148.0 
{E*=31.5) MeV to £;cm=176.0 (£;*=61.5) MeV (see the 
histograms). At high incident energies, the mass distri- 
bution has a Gaussian-like shape centered at the sym- 
metric mass division, whereas the mass-asymmetric fis- 
sion fragments dominate at low incident energies. The 
mass-asymmetric fission produces nuclei in the vicinity 
of the doubly-closed shell nuclei, ^^^Pb and '^^Ni. The 
trend of the experimental data, i.e., the incident energy 
dependence of MDFF, is well reproduced by the calcula- 
tion. The mass-asymmetry with ^h = 200 at sub-barrier 
energies is also well reproduced. 

In Fig. [2] we also plot the fusion-fission events by the 
filled histograms. Apparently, the compound-nucleus fis- 
sion has a mass-symmetric shape, and the observed mass- 
asymmetric fission dominated at the low incident energies 
is classified as QF. The strong energy dependence of the 
MDFF can be understood in terms of the orientation ef- 
fect on the fusion and QF. The collision on the polar side 
have a large probability to disintegrate as QF, whereas 
the collision on the equatorial side have a larger fusion 
probability. The standard deviation of the spectrum for 
FF is cr,„=37.2 u at i^c.m.^lTG.O MeV. This value is far 
smaller than that of the experimental spectrum, which 
also shows the Gaussian-like shape, indicating that there 
is a significant contribution of QF even at this low en- 
ergy. The calculation also suggests that the measured 
mass-symmetric fission fragment has another origin than 
the compound nucleus fission. Such an event is defined 
as a deep quasi-fission process (DQF) as discussed in Ref. 

The results for the '^'*S -I- ^ssy reaction are shown in 
Figs. [3] and S) In addition to the effects of the static 
deformation of 238^ g^j^^j ^j^g couplings to the 3~ state 
in ^■^^U, we also take into account the 2+ state at 2.13 
MeV in ^''S in the CC calculation. Figure H shows 
the MDFF for five incident energies between £^c.m=148.0 
(£;*=36.2) MeV and Sc.m^lVO.O (£;*=58.2) MeV. The 
Bass barrier for this reaction is Veass = 161.1 MeV. The 
meaning of each curve in Figs. [3] and [3] is the same as 
in Figs. [1] and [2l The energy dependence of MDFF is 
qualitatively the same as in the "^^S + ^ssy reaction, and 
the calculations again reproduce well the experimental 
data. 



B. Reaction of 3081+^^* 



U 



the kinetic energy has —z direction, while in the previous 
study, we started the Langevin calculation from a suffi- 
ciently large distance between the target and projectile 
and the orientation dependence was not considered. The 
previous results of smaller fusion probabilities are due to 
the loss of the kinetic energy in the approaching process 
by the friction. 



We next analyze the '^"Si-I-^'^^U reaction, for which the 
MDFF and the fission cross sections have been measured 
by the JAEA group [13|. The present calculation for 
the capture cross section (Tcap, fusion cross section Cfus 
and the mass-symmetric fission cross section aA/2±20 ^-re 
shown in Fig. [5] by the dashed-dot, solid and dashed 
curves, respectively. For the capture cross sections, sim- 
ilar parameters for CC are used as in the ■^^^'^^S -I- ^38^ 
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FIG. 3. (Color online) Same as Fig. 3, but for the ^''S+^^^U 
reaction. The experimental data are taken from Ref. [l7| . 



reactions, except that the couplings to the vibrational 
state in ^°Si are not taken into account. The present cal- 
culation for the capture cross sections reproduces quite 
well the fission cross sections down to the sub-barrier re- 
gion below the Bass barrier, Veass = 141.1 MeV, better 
than the previous model calculation [1^. The predicted 
fusion cross sections fjfus are about 5 to 7 times larger 
than the previous model 19] due to the same reason as 
in the '^^'•^^S-f^'^^U reactions discussed in the previous 
subsection. 

The calculated results for MDFF are shown in Fig. [S] 
by the histograms for eight incident energies i?c.m (E*) 
between 129.0 (35.5) MeV and 169.0 (75.5) MeV. We 
use the same parameters for the trajectory calculation as 
those for the 36,34g _|_ 238u reactions. A significant dif- 
ference in the measured MDFF between the ■^°Si -I- ^ssy 
and the 34,36g _|_ 238y reactions is the mass-asymmetry 
in QF at subbarrier energies. The mass-asymmetry for 
the 3°Si -I- ^^^U system is ^hMl ~ 178/90 as shown in 
Fig. M whereas AhMl is 204/68 and 200/74 for the ^48 
-I- 238^ g^j^^j ^Yie ^^S + 238u reactions, respectively (see 
Figs. Hand 15). The QF fragments Ah/Al = 178/90 do 
not fit any shell closure of neutron rich nuclei. A pro- 
duction of these nuclei is therefore not associated with 
the local minimum of the potential energy attained by a 
larger binding energy of the nascent fission fragments at 
the scission point, but should originate from dynamical 
aspects in heavy-ion induced fission. This feature is well 
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FIG. 4. (Color onhne) Same as Fig. 4, but for the ^''S-H^'^^U 
reaction. The experimental data are taken from Ref. [17l |. 



reproduced by the present calculation, in the entire en- 
ergy range from above-barrier to sub-barrier. The stan- 
dard deviation of the measured MDFF decreases from 
the highest incident energy down to the low energy of 
Ec.m. = 139.4 MeV, then the value suddenly increases at 
the sub-barrier energies of 134.0 and 129.0 MeV. These 
trends are also reproduced by the calculation. 

The calculated FF events are shown by the filled his- 
tograms in Fig. [S] It has a mass-symmetric distribution. 
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FIG 



5. (Color online) Same as Fig. 3, but for the 30si_^238 



U 



reaction. The experimental data are taken from Ref. [IS 



and the standard deviation does not change through all 
the energy points. Therefore, the sudden broadening of 
the measured MDFF at sub-barrier energies can be at- 
tributed to the enhancement of the mass-asymmetric QF 
at ^hMl = 178/90. 



IV. ANALYSIS OF THE STRUCTURE IN THE 

MASS DISTRIBUTION OF FISSION 

FRAGMENTS 

In the previous section, we have shown that the new 
model nicely reproduces the energy as well as the sys- 
tem dependence of MDFF. In this section, we discuss 
the origin for the difference in the mass-asymmetry in 
QF between the ^"Si + ^^^U and ^^S + ^^^^U systems 
based on our model. 
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FIG. 6. (Color online) Same as Fig. 4, but for the ^'''Si-h^^^U 
reaction. The experimental data are taken from Ref. [la |. 



A. Potential energy surface along the scission line 

Let us first consider the potential energy surface for 
the two systems. Generally, the shape of MDFF from 
a compound nucleus is affected by the landscape of the 
potential energy surface, especially near the fission sad- 
dle point and the scission point [44, 45]. On the other 
hand, for QF, the saddle point of the composite system is 
not passed by the system during the evolution of nuclear 



shape. A plausible explanation for the mass-asymmetry 
in QF may have something to do with the potential en- 
ergy at the scission point as the exit point of the reaction 
process. Therefore, we focus on the landscape of the po- 
tential energy surface near the scission line. 

Figures [7ja) and (b) show the adiabatic potential en- 
ergy surface near the scission point for the nuclei ^^*Sg 
(30gj _,_ 238u) and 274hs (36s + '^^^V), respectively. The 
dashed and solid curves correspond to the potential Vld 
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FIG. 7. (Color online) Adiabatic potential energy surfaces 
near the scission point for the nuclei ^®*Sg (a) and '^^''Hs (b). 
Vld and Vld + i^shcii with 5 = 0.22 are represented by the 
dashed and solid lines, respectively. The red and blue lines 
denote the potential with e — 0.35(z = 2.35, 5 — 0.22) and 
e — 1.0{z = 1.5, S = 0.22), respectively. 



and VhD + ^fhou '^itl^ deformation oi S = 0.22, respec- 
tively. This value of the deformation parameter S is cho- 
sen since it reproduces well the available fission data as 
discussed in Refs. [4^|43]- As we mentioned in Sec. II 
B, the value of e evolves from e = 1 towards e = 0.35 as 
a system approaches the compound-nucleus shape and 
fissions. The blue and the red lines show the potential 
surfaces for the configurations corresponding to these val- 
ues of e (e=l and 0.35), that is, (z, 5) = {1.5, 0.22) and 
(z, S)—{2.35, 0.22), respectively. The sohd and dashed- 
dot curves in Fig. [5] represent the corresponding nuclear 
shapes with a mass symmetry a =0. Notice that for a 
faster process like QF, the system does not have enough 
time to reach e — 0.35 and e is maintained about 1.0 still 
at the scission point. 

Even though the shell energy E^y^^^^ has a strong local 
minimum at A ^ 208 (which corresponds to ^°*Pb) ow- 
ing to the strong shell structure, the local minimum of 
Vld + -E'shoii appears at A ^^ 186 due to the parabolic 
shape of Vld around A - 135 for both ^esgg and ^'^'^Hs. 
One also finds another minimum at A ~ 85 for both the 
nuclei. If the mass division were determined at the scis- 
sion point, the mass- asymmetric QF should thus have 
appeared at A ~ 86 and 188 for the reaction ^^S + ^^^U 
and at A -- 84 and 184 for the reaction ^"Si -I- ^^^U. 
These are not realized in the measured MDFF, however. 
The two peaks observed in the experiment are located at 
approximately A ~ 74 and 200 in the reaction of ■^^S -I- 
^^^U. Evidently, a prediction based only on the poten- 
tial energy surface at the scission line is insufficient to 
explain the mass asymmetry in QF, indicating that the 
dynamical process plays an essential role to describe QF. 
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FIG. 8. (Color online) The nuclear shapes near the scission 
point. The dashed-dot and solid lines represent the nuclear 
shapes corresponding to the two curves in Fig. 9, that is, 
e = 0.35 {z = 2.35, 5 = 0.22) and e = 1.0 {z = 1.5, 5 = 0.22), 
respectively. 



B. Analysis of the reaction dynamics using 
probability distribution 

We next discuss the probability distribution of the sys- 
tem in the deformation space. To this end, we segment 
the coordinate space with Az = 0.10, A(5 = 0.06 and 
Aa = 0.06. We define the distribution as an ensemble 
of Langevin trajectories. That is, we follow a trajectory 
as a function of time, and we increase the event number 
at each segment when the trajectory passes through that 
segment. By generating many trajectories, we construct 
a distribution of events on the deformation space. 

Figure |9] shows the distribution so constructed on the 
z-A plane for the reaction ^osi + ^ssy at E* = 35.5 MeV 
with £ = 0, and = Q. We consider the time evolution of 
the nuclear shapes as well as the mass number only for 
one of the fragments, that is, the one corresponding to the 
target-like nucleus. Even though this distribution is not 
a probability distribution, because it is not normalized 
in the whole space, we call the distribution "probability 
distribution" . From such probability distribution shown 
in Fig. [SI we can understand the overall trend of the 
dynamical evolution of nuclear shape, where the vertical 
axis of the figure means the number of events entering in 
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FIG. 9. (Color online) The probability distribution con- 
structed as an ensemble of Langevin trajectories on the z — A 
plane for the reaction ^°Si + ^^®U at E* = 35.5 MeV with 

£ = 0, e = 0. 



FIG. 10. (Color onhne) Same as Fig. 14, but for for the ^''S 
+ ^^*U reaction at E* = 40.5 MeV with ^ = 0, 61 = 0. 



(a) "Si+^'U, E*= 35.5 MeV 



the z — A plane. 

Two remarkable peaks are visible in Fig. |9l The peak 
located around [z^A) ~ (1.3,238) corresponds to the 
touching point. All the trajectories start at this point, so 
that the number of events at this segment is large. The 
broadly distributed events in the smaller- z region corre- 
spond to the trajectories trapped in the large-(5 region 
as seen in Figs.?? (a) and (c). They are the trajectories 
that enter into a compact region defined by z < 0.5 and 
75 < A < 150, which correspond to afus in Fig. [5] We 
also see the yields which leads to scission at ^ ~ 178 and 
135 dX z — 2.5, which correspond to the peaks of MDFF 
shown in Fig. [6l The peak Sit A ^ 178 originates from 
QF, while the peak at A ~ 135 originates largely from 
the FF, but also involves the component from DQF. 

The probability distribution for the ^^S -I- ^^^U reac- 
tion at E* = 39.5 MeV is shown in FiglTUl The spectrum 
is dominated by the mass-asymmetric QF. Similarly to 
the ^°Si-|-^'^®U reaction, the large peak corresponds to 
the touching point. Almost all the trajectories move 
quickly to the direction of binary decay as QF, and those 
which approach the compact nuclear shape with small z 
value, corresponding to Ufus in Fig. [l] are significantly 
diminished. QF generates fragments around A ^ 200 
(at z = 2.5) in this calculation, explaining the observed 
mass asymmetry of MDFF in Fig[2]at the lowest incident 
energy. 

The probability distributions projected onto the z — 5 
plane are shown in Fig. [TlTa) and (b) for the reac- 



tions of -^"Si 



238 



U and 36 s 



238 



U, respectively. The 



dashed curves denote the ridge lines. The positions at 
the strongest yield correspond to the touching points. 
In Fig. [TIT a) , the distribution crosses the ridge line and 
moves to the large-(5 region. It spreads in the small-z re- 
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(b) "S+'^'U, E'= 40.5 MeV 
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FIG. 11. (Color online) The contour maps for the probability 
distributions on the the z — 5 plan, for the reactions of '^"Si 
+ 2^*U at E* = 35.5 MeV (a) and ^"^S + ^^*U at E* = 40.5 
MeV (b) with I = Q, 6 = Q. The dashed curves correspond to 
the ridge lines. 



gion and 5 value distribute widely from —0.1 to 0.7. Af- 
ter surmounting the ridge line and entering in the com- 
pact nuclear shape region, a mono-nucleus with small 
z-shape is formed, but its deformation fiuctuates signif- 
icantly due to the thermal fiuctuation. The distribution 
for z < 0.5 and 5 < 0.2 corresponds to the formation of 
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the compound nucleus. The ^^S + ^38^ reaction shown in 
Fig- HUb) has a much smaller chance to form the mono- 
nucleus as compared to the reaction ^°Si + ^'^^U, because 
of the smaller probability to overcome the the ridge line. 
Instead, the system quickly disintegrates as QF by form- 
ing fission fragments. 

We can clearly observe that the fission fragments for 
the reaction ^°Si -I- ^^^U are deformed with -0.2 < 6 < 
0.5, because they include various events originating from 
QF, DQF and FF. On the other hand, for the ^^S + ^^^V 
reaction, the value of the deformation parameter of the 
fission fragments are restricted mainly to —0.1 < S < 0.2, 
primarily because QF is the dominant channel. 

From the analysis of the probability distribution, it is 
clear that the ^^Si -I- ^^^U and ^^^S + ^^^U reactions at 
low incident energies are governed by significantly differ- 
ent reaction processes, that are reflected in the peaks of 
MDFF. For the former reaction, the peak at A '^ 178 (the 
asymmetric peak) originates from QF, while the peak at 
A ~ 135 (the symmetric peak) originates both from DQF 
and FF. On the other hand, the mass-asymmetric peak 
for the latter reaction is located at A ^ 200, originating 
predominantly from QF. 



V. SUMMARY 

We developed a new dynamical model to describe 
heavy-ion induced fission, in which the effects of static 
nuclear deformation of a target nucleus are taken into 
account by considering all the orientation angles of the 
symmetry axis of the target nucleus. The orientation ef- 
fects are included both in the barrier penetration process 
and in the evolution of the nuclear shape. The former 
process is described with the coupled-channels model. 
After the nuclear contact point, we switch to the dynam- 
ical calculation starting at the touching point assuming a 
nose-to-nose configuration. The angle dependent touch- 
ing distance ztouch is introduced in order to effectively 
simulate the orientation effects in the evolution phase 
of the nuclear shape. With this model, the calculation 
could be extended to energies below the Coulomb barrier 
for the first time. 

In spite of the simplified assumptions, the calcula- 
tion reproduced the measured MDFF for the reactions 
of 36,34g_(_238u ^^^ 30gj_,_238u^ where a large variation of 
the distribution with respect to the incident beam energy 



are experimentally observed. By analyzing the Langevin 
trajectories, we could distinguish three different fission 
processes, QF, DQF and FF. The mass-asymmetric fis- 
sion in 368-^238^ ^^ Al/Ah = 74/200 observed at low 
incident energies are from QF. The mass- asymmetric fis- 
sion with Al/Ah = 90/178 in 30si-|-238u observed at 
sub-barrier energy is also from QF, whereas the peak at 
the symmetric fission of yl=135 indicates that the fission 
occurs with some time delay and originates from FF and 
DQF. The calculation suggests the formation of a mono- 
nucleus with a small z value but with a large collective 
fiuctuation on the deformation axis of 5. A relatively 
long-life mono-nucleus is formed by the potential pocket 
appearing inside the ridge line. During the process of 
surmounting the ridge line, the mass asymmetry (a) of 
the system moves to the symmetric region. This can be 
an account for the measured mass asymmetry in QF in 
30gj_,_238u^ The trajectories for ^'^S-I-^^^U can not cross 
the ridge line and are directed to QF without forming 
the mono- nucleus, thus the mass-asymmetry is close to 
that of the entrance channel. 

In this model, one can determine the fusion probability 
by selecting trajectories which enter the fusion box. The 
reproduction of the experimental MDFF in this model 
can be the ground to support the calculated fusion prob- 
ability. Furthermore, the generalized formula proposed in 
this model has a potential to simulate any kind of heavy- 
ion induced reactions in the approaching phase, such as 
a nucleon-transfer reaction, and to predict cross sections 
for the production of new nuclei. 
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